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The proton form factors at large momentum transfer are dominated by two contributions which are 
associated with the hard and soft rescattering respectively. Motivated by a very active experimental 
form factor program at intermediate values of momentum transfers, Q 2 ~ 5 — 15 GeV 2 , where an 
understanding in terms of only a hard rescattering mechanism cannot yet be expected, we investigate 
in this work the soft rescattering contribution using soft collinear effective theory (SCET). Within 
such a description, the form factor is characterized, besides the hard scale Q 2 , by a hard-collinear 
scale QA, which arises due to the presence of soft spectators, with virtuality A 2 (A ~ 0.5 GeV), such 
that Q 2 3> QA 3> A 2 . We show that in this case a two-step factorization can be successfully carried 
out using the SCET approach. In a first step (SCET/), we perform the leading-order matching of 
the QCD electromagnetic current onto the relevant SCET/ operators and perform a resummation 
of large logarithms using renormalization group equations. We then discuss the further matching 
onto a SCET// framework, and propose the factorization formula (accurate to leading logarithmic 
approximation) for the Dirac form factor, accounting for both hard and soft contributions. We also 
present a qualitative discussion of the phenomenological consequences of this new framework. 



The study of the nucleoli form factors (FFs) is one of the central topics in hadronic physics ( for 
recent reviews see, e.g. , Refs. [IH3])- Substantial progress has been achieved in this field over the past 
decade, mainly thanks to new experimental methods, using polarization observables, which allow for 
precise measurements of the FFs. The results for the proton FFs, obtained over the past few years at 
JLab [4-7 up to a momentum transfer Q 2 ~ 8.5 GeV 2 , considerably boosted our knowledge about the 
distribution of the electric charge inside the proton. A substantial program to extend the measurements 
of the nucleon FFs up to Q 2 ~ 15 GeV 2 in the spacelike region will be performed in the near future at 
the JLab 12 GeV upgrade. In parallel, the PANDA Collaboration at GSI is planning to carry out precise 
measurements of the proton FFs at large timelike momentum transfers, up to around 20 GeV 2 , using the 
annihilation process p + p — > e + + e~ [8 . These experiments will provide us with precious information 
on the FF behaviors in the region of large momentum transfers. 

On the theory side, an understanding of the nucleon FFs at large momentum transfers, both spacelike 
and timelike, from the underlying QCD dynamics, still remains a challenge. At present, the FF behavior 
for moderate and large values of Q 2 is still not well understood and an adequate description, allowing for 
quantitative predictions, is absent. 

The leading power behavior of the FFs was studied a long time ago using the QCD factorization 
approach ( see, e.g., [21 HO] and references therein). It was established that the dominant contribution 
can be represented by a reduced diagram as shown in FigjT] In this figure, the hard blob describes 
the hard scattering of quarks and gluons with virtualities of order Q 2 . Such a hard subprocess can be 
systematically computed in perturbative QCD (pQCD) order-by-order. The soft blobs, denoted by 
describe the soft, nonperturbative subprocesses, and can be parametrized in terms of universal matrix 
elements known as distribution amplitudes (DAs). Such a picture suggests the well known factorization 
formula for the Dirac FF 



I. INTRODUCTION 
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FIG. 1: Reduced diagram describing the hard scattering picture 



and predicts the scaling behavior 

Fi~^xpnQ/Af , (2) 

where the logarithmical corrections can be systematically computed order-by-order. 

Unfortunately for the Pauli FF F 2 this approach cannot provide such a systematic picture and suggests 
only the power estimate 



A> 



*i~o«r (3) 



Almost simultaneously, it was found that the picture described by Fig[T]is not complete. In Ref. [TT] 
it was demonstrated that the exchange of soft quarks between initial and final states may also produce 
contribution of order 1/Q 4 times logarithms. In Refs. [321 [13] all such contributions were computed with 
the leading logarithmic accuracy at 2 and 3 loops. Using these results it was assumed [13] that these 
"nonrenormalization" logarithms probably can be resummed to all orders into an exponent similar to 
the well known Sudakov logarithms [14] . However, this effect was ignored in many later publications. In 
particular, in Ref. |10j it was suggested that such contributions could be strongly suppressed due to those 
Sudakov logarithms and therefore can be ignored at large values of Q 2 . 

At the same time, many phenomenological studies of the hard rescattering picture support the conclu- 
sion that in the region of moderate Q 2 ~ 5 — 10 GeV 2 the factorization approach expressed by Eq. 
cannot describe the data properly ( see, for instance, [H]). Moreover, existing data for the FF ratio 
F2/F1 measured up to Q 2 = 8.5 GeV 2 [7] also do not support the expectation of Eq. ^ which assumes 
that in the asymptotic region Q 2 F2/F\ ~ const. Therefore, it was suggested that the so-called Feynman 
mechanism |16) , associated with the scattering of the hard virtual photon off one active quark, dominates 
the nucleon FFs at moderate values of Q 2 . The other spectators remain soft and therefore very often 
such scattering is associated with the soft overlap of the nucleon wave functions. 

Such a picture is supported by different phenomenological approaches, such as QCD-motivated models 
for the hadronic wave functions [T7HT9] , QCD sum rules [20l [2Tj and light-cone sum rules [22| [23] . The 
Sudakov suppression in this case is always assumed to be relatively small. The aim of the present work 
is to develop a systematic approach for the specific soft contribution described first in Ref. [TT], and 
formulate it through a factorization theorem. We apply the effective theory approach, known as soft 
collinear effective theory (SCET), in order to describe contributions from different regions of virtualities 
in the diagrams. 

The effective theory is a very convenient tool in this case because soft rescattering is characterized 
by subprocesses which exhibit different scales : a hard rescattering involving particles with momenta 
of order Q 2 , hard-collinear scattering processes with virtualities of order AQ, and soft nonperturbative 
modes with momenta of order A 2 . Therefore one has to perform a two-step matching procedure in order 
to perform full factorization of such a process. 

Following this scheme we obtain that the full description of large-Q 2 behavior of the nucleon FF Fi is 
given by the sum of two contributions associated with the soft and hard rescattering picture: 



F^+Ff). (4) 

The hard rescattering part F^ is well known and described by fl. One can expect that the soft contri- 
bution can also be presented in a factorized form but with the more complicated structure reflecting the 
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presence of different scales. Performing the leading logarithmic analysis of the leading power contribution 
(~ 1/Q 4 ) we demonstrate in this work that the corresponding soft term can be presented in the following 
form, 

//>oo /> />oo 

Dy^ivi) J du>idu> 2 J'faMQ) J D Xi V(Xi) J dv x dv 2 J(x is ViQ)S{u) U (5) 

which can be interpreted in terms of a reduced diagram as in Fig[2] This result involves a hard coefficient 
function H, and two hard-collinear jet functions J and J' which can be computed in pQCD. They describe 
the subprocesses with hard momenta and hard collinear momenta respectively. The nonperturbative 
functions and S describe the scattering of collinear and soft modes. The convolution integrals in 
Eq. ([5| are performed with respect to the collinear fractions Xi and yi, and with respect to the soft 
spectator fractions Wj, ~ A. 




FIG. 2: Interpretation of the soft rescattering as a reduced diagram 

In the case of the Pauli FF F 2l we can also perform a factorization of the soft-overlap contribution but 
only partially, separating the hard modes with momenta of order Q 2 . The full factorization is problematic 
due to overlapping integration regions corresponding with soft and collinear contributions, which lead 
to well known end-point singularities in the convolution integrals. However, such a partial result can 
be used to carry out a phenomenological analysis of the FFs in the region of intermediate Q 2 values. 
Such a region corresponds to momentum transfers where Q 2 is large enough, allowing us to perform a 
power expansion, but where the second, hard-collinear scale ~ AQ is still relatively small, so that one 
expects the dominance of the leading power asymptotic term. Such a situation may indeed be relevant 
to interpret existing data and planned experiments. 

The specific feature of the factorization for the soft-overlap contribution is the presence of the Sudakov 
logarithms which can be ressummed using the renormalization group in effective field theory. It was 
suggested, see e.g. [24] , that these logarithms could play an important role in the timelike region 
providing an enhancement of the timelike FFs compared to the spacelike region (the so-called K factor) . 
Within the factorization picture such an enhancement can be clearly studied in a model independent way. 

Our paper is organized as follows. In Sec. [TT] , we consider as an example the analysis of the dominant 
regions for certain Feynman diagrams and demonstrate the existence of the soft spectator contribution 



at leading power (in the hard scale Q) for both the Dirac and Pauli FFs. In Sec. Ill we discuss the 



factorization scheme for such contributions, perform the leading-order matching between QCD and the 



SCET, and perform a resummation of large logarithms. In Sec. IV we discuss the SCET power counting 
in 1/Q and derive the factorization formula In Sec. |vj we perform a first qualitative discussion of 
the phenomenological consequences following from our results. In Sec. |VI| wc summarize our findings. 



II. SOFT RESCATTERING MECHANISM: EXAMPLES 



In this section we consider specific examples of soft rescattering contributions. For the Dirac FF our 
analysis overlaps with results of the work of Ref. [TT] , whereas for the Pauli FF this is discussed here for 
the first time. 
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In our consideration we use the Breit frame 

q = p' -P = Q{^- , n = (1,0,0,-1), n= (1,0,0,1), (nn) = 2 
and define the external momenta as 

1 



P=Q 



n m N n 
2 + ~Q2' 



p=Q 2 + ^Tr Q = Q 2 



i + Wi 



Am 2 N 



Q 2 



= Q + 0(m N /Q ), 



(6) 



(7) 



2(pp>) = Q 2 + ^*Q 2 , 
where m N is the nucleon mass. For the incoming and outgoing collinear quarks we always imply 



Pi = XiQ- +pj_ t 
with the transverse momenta 



Q 



Vi 



m" N \ n 
Q J 2' 



(8) 



(9) 



pI-p'I-a 2 , 

and where Xi and x\ denote fractions of the corresponding momentum component. In what follow we 
shall use the convenient notation Xi = 1 — x^. We also use the following notation for scalar products 



and Dirac contractions 



(a ■ n) = a + , (a • n) = a_ . 



£V7 M =1> = P- 



(10) 



(11) 



Nucleon FFs are defined as the matrix elements of the electromagnetic (e.m.) current between the 
nucleon states: 



<pV£ m .(o) \p) = n(p') 



2m n 



N(p), 



(12) 



with nucleon spinors normalized as NN — 2rriN- We also use a standard normalization for particle states: 

{p',s'\p,s) = (2ir) s 2E8 ss ,6(p-p'). (13) 

In what follows, we shall compute the Feynman diagrams which provide contributions to the nucleon 
FFs. The component of interest for our calculations is the soft matrix element describing the overlap of 
the partonic configurations with the hadron state. In the case of the FFs such overlap is described by 
DAs. In the case of the nucleon, the corresponding leading twist DAs can be defined as 



4<0| MQ ^[A 1 n] U(9 VF J [A 2 n]d CT ^ fc [A 3 n]|p) = ^- J Dx t e - ip +^ XtXi ^{xi) 



(14) 



where 



q a W [x] = q a (x)P cxp l^ig J dt (n ■ A) (x + tn) | , (15) 
and the measure reads Dxi — dxidx2dx 3 5(l—xi—x 2 —x 3 ). The function &(xi) can be further decomposed 



as 



= V{ Xi ) P+ [a# C] af3 [j 5 N+] a + A(x l ) P+ [^l 5 C} a/3 [N+] c 



(16) 
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The large component N + of the nucleon spinor is defined as 



4 ' 



(17) 



and C is the charge conjugate matrix C : C 1 j ll C = — 7^ . The nucleon DA ^f(xi) is shown by the 



"^C = - 7 J. 

soft blobs in FigTlJ For simplicity we restrict our consideration to the proton state. In what follows we 
always assume that in pQCD diagrams the first and second top lines correspond to u quarks. Assuming 
projections on the leading twist DAs, we can considerably simplify certain considerations substituting 
instead of DAs on-shell quark spinors. Such substitution is possible because at leading-order (LO) power 
accuracy we can neglect the small components in the external quark momenta: 



XiP, p 



ViP 



(18) 



and assume that external collinear quarks are on-shell. This is possible because the leading twist projec- 
tors [T x C], [T x N + ] a in |i~4|) satisfy to following relations 



P = pT x =pT d x N+ ~0, 



(19) 



that are compatible with the free equation of motion for quark spinors and allow us to use on-shell 
spinors in the intermediate calculations. The contribution to the physical amplitude can be obtained by 
resubstitution of the quark spinors by the hadronic matrix element (14). 



A. Soft rescattering contribution for the Dirac FF F\ 

Consider, following Ref. [IT], the diagram in FigH The incoming and outgoing particles must have 
invariant mass ~ A 2 in order to overlap with nucleon states. This is guaranteed by the momenta in 
Eq. Q. The interactions between the external quarks are soft and described by DAs. For simplicity, the 
corresponding soft blobs are not shown in Fig. [3j One can easily find that 



3 k i i 


'1 p, 1 




= 

g 

p 2u p k 2 a p' 2u = 


Pi P ' 


c 

P~> ° 
c 




P 4g P 



3d 



P3 



FIG. 3: The simplest diagram with soft exchanges. 



Pig =Pl - fa, Pig =P~PS - h - k 2 , p 3 g = fa -P'%, Pig = fa + fa ~ p' + P3, (20) 

P2u=p-P3-fa, P3d =P~ fa - k 2 , p 2u =p' -P 3 - fa, Ptid = V' - fa - fa- ( 21 ) 

The analytical expression for the diagram of Fig. |3j where the quark line with momenta p 3 and p' 3 
represents a d quark, reads 



Df, = C 



cf k 1 d 4 k 2 



1 d(p' 3 ) l a (P' -fa- fa) r (p-fa- fa) l p d(p 3 ) 

[k 2 2 - m 2 ] [k\ - m 2 ] (p/ _ kl - k 2 f (p-fa- k 2 f (pi - fa) 2 {kt - pi) 2 

l (Pih l (fa + m ) l Ju (Pi) ^{p'2) 1 1 {p' -P 3 ~ fah a (fa +m)~fP(p-p 3 - fa)^ u(p 2 ) 

(P-Pi- faf (P' ~P 3 - faf (fa + fa ~P' +P' 3 f (P~P3-fa- faf 



(22) 
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where the numerical factor C accumulates all color factors and vertex and propagator factors, and m 
denotes the quark mass 1 . We write on-shell quark spinors instead of projectors on the nucleon DA as it 
was described above. 

According to the factorization expressed by Eq. , one could expect that dominant integration regions 
(providing contributions of order A 4 /Q 4 ) can be described as follows: 

hard region: fcf - Q, k 2 - Q 2 , (23) 
collinear-p region: k { : (kn) — Q, (kn) — A 2 /Q, k± — A , kf — A 2 . (24) 

collinear-p' region: k, : (kn) ~ Q, (kn) ~ A 2 /Q, k x ~ A , k 2 ~ A 2 , (25) 

Then factorization formula ([!]) implies that the general structure of any 2-loop diagram can be interpreted 
as 

D = * * T (2) * * + * T (1) * * + * * T (1) * 

+ *( 2 ) *T (0) * * + * *T (0) * * (2) , (26) 

where 

*W=VW**, = V (1) * V (1) * *, (27) 

denotes the convolution of the collinear evolution kernel of order i with DA \&. Such contributions 



related with the collinear regions (24 1 and (25 1. The hard kernels T^ ' 1,2 ' denote the contributions 
to the hard coefficient function in LO, next-to-leading order and next-to-next-to-leading order order, 
respectively. 

However, this description is not the full answer at the leading-order accuracy in 1/Q. There is one 
more region, which cannot be interpreted in the form of the reduced diagram in Fig. [I] and is defined as 
the soft region: 

fcf~A, k 2 - A 2 . (28) 



Let us compare the power contribution from this region with the contribution from the hard region ([23j) . 
The latter provides 

Dfl - Q 8 ^ ~ Q'j^o £r 2 6 g/j, ~ ^ ^ gr 3 & (29) 

where 1^ denote some scaleless Dirac structures, and the factor ~ Q s arises from the measure. The term 



with 7j* in ( 29 1 reflects the requirements of one transverse index ~ 7j_. In order to arrive at the formula 



(29) we also used the decomposition of the quark spinors onto large and small components: 



m = s + fft, £ = m", % = m^f, (30) 

qXPi) = & + Vi, 6 = -^<l(Pi)i Vi = -j-QKPi)- ( 31 ) 
One can easily obtain that the small component n is suppressed relative to the large component £ as 

V ~ C/Q- (32) 



1 For simplicity we do not show explicitly the color indices. The quarks mass in written only in the propagators where it 
can be relevant. 
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Consider now the contribution from the soft region. Neglecting in the denominator (22 1 by small terms 
one obtains: 



Den : 



[kl - ™ 2 ] [*? - ™ 2 ] {yiy'i [-V ■ (fci + fc 2 )] 2 [-2(fe! -p')] 2 } 

[ Xl x\ [-2p ■ (h + k 2 )f [-2 (A* • p)] 2 } ~ A 2 A 2 (AQ) 8 , 



where recall, Xi = 1 — Xj. Therefore 



„t*\ . a Num , . 

£)W ^ A 8 A s 

Den 



Num 



A 2 A 2 (AQ) 



(33) 



(34) 



From Eq.(34) we see that the numerator must contain the soft scale at least in the power ~ A 4 or higher. 

us in tl 
ne we o 

2) Y l (p-k- k) i p d( P3 ) ~ d( P ' 3 ) rp'i^ d(p 3 ) 



We next compute the largest terms in the numerator. Neglecting the small momenta and the small 
spinor components in the d-quark line we obtain: 



- ki - k 2 



Then the second u-quark line can be rewritten as 

p"V u(p' 2 ) f ip' -p' 3 - fci)7 a (fc 2 +m) 7 ^(p-p 3 - A;i)7 J u(p 2 ) ^ € 2 l i fap'(h + 'm)pkij : > £ 2 
The product of the w-quark lines yields 

CW 4 (fci + m) l ] W 2 l%p'{k 2 + m)phy 6 

~2(p' • fci)2(p • fci) CW l (fci + m) j j h C 2 f(k 2 + m) 7 j 6- 

Therefore we obtain 

Num = 2(p' ■ fcx)2(p ■ fci) |i 7 « (fc! + m) 7^1 &y'(*2 + m)V 6 & 

~ Q 2 A 4 gr^ gr a & fh^a- 



(35) 
(36) 



Substituting this into ( 34 ) yields 



Q 2 A 4 



DW~A 8 , 

Pi A 2 A 2 (AQ) 



gr a & ^7^3 



(37) 



(38) 



(39) 



i.e. we obtain the same power of Q as for the hard region in (29). Therefore we established that the soft 
region is the additional relevant region which is not accounted for in the factorization formula 0. In 
Ref. [12] all diagrams with the soft spectator quarks have been computed with the leading logarithmic 
accuracy. Their sum does not cancel providing some nontrivial answer. Hence we can avoid consideration 
of such a possibility. 

Consider now the whole expression for the soft contribution: 



d(;1 = AC / d%d% [<# ] a 



303 



(fci+m) (k 2 +m) a2 p 2 



[14 



, 2 ] [k\ 



yi yl[-2p' ■{k l + k 2 )f[-2{k 1 - P ')\ 



7^U [7 J 6] & M 



ft 



Xl x\ {-2p-(k 1 +k 2 )Y[-2{k 1 -p)] 



(40) 



Each expression in the square brackets describes some subprocess involving the particles with appropriate 
virtualities and momenta. We consider them term by term. The factor 



(ki + m) (k 2 + m) a2 p 2 

7 7 °ift 

[fef — m 2 ] [kl — m 2 ] 



(41) 
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describes the propagation of the soft spectator quarks and includes only soft particles with kf ~ A 2 . This 
term can be associated with the soft part of the diagram. The factor 



yi yl[-2p' ■(k 1 +k 2 )] 2 [-2{k 1 -p')\ 



(42) 



describes the transition of two soft spectator quarks and one active quark into three collinear quarks. 
It is described by the subdiagram with the two-gluon exchange. As one can see, all propagators have 
virtualities of order fc, • p' ~ QA and all involved momenta have a large component ~ Q along the p' 
direction. 

In a the similar way one can describe the second subprocess given by 



Xl xl \-2p-{k 1 + k 2 )f[-2(k 1 -p)] 



(43) 



The difference from the previous case is only in the involved momenta. They have large components ~ Q 
along the p direction. 

The simple vertex factor [7^ ] a3 g 3 can be associated with the hard scattering vertex of the subprocess 
7* + d(p) —> d(p'). It is clear that this subprocess in general involves particles with large momenta of 
order Q 2 . 

Taking into account the different virtualities of the particles: A 2 -C QA <C Q 2 one can try to factorize 
the whole result of Eq. (40) in accordance with the described subprocesses. In order to do this we 
introduce the Sudakov decomposition 



{h ■ n) 



kf, (ki ■ n) = k t , dki = -dkfdk i dk^ 



(44) 



and rewrite Eq. ( 40 ) as 



Ml La, / dkfdk 



C / dki 



1 W] ai W] a2 [&U 

Q 3 vivl (kt + k+y [-k+] 



[ k\ + m J ( k 2 + m ) 

V /oift V /a 2 /5 

[k\ — to 2 ] [k 2 — m 2 ] 

Q 3 xix\ (fcf + k 2 ) 2 [-/cf] 



(45) 



This equation almost represents the required form. To make it more obvious Eq. ( 45 ) can be rewritten 
as: 



where we introduced 

S( 0J 9) {u^Vi) = C J dk^dk 2 dki 2 ±S(u> 1 - kf)S(uj 2 - k%)8{v\ - k^)5{v 2 - k 2 ) 
(ki+m) lk 2 +m) 



(46) 



CC202 



[kl 



The two functions J[ > and Jtp\, which we will refer to as jet functions, read 



1 1 



1 



? 3 2/12/1 ( W l + W 2) 2 [-Wl 



[3]« 3 [^l] Q2 [Ci7l] Ql 



716 fl 7i^i 

J P2 L 



(47) 

(48) 
(49) 
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where the index in brackets denotes multi-index: (a) = {ai,a2,a 3 }. These functions describe the 
scattering of the particles with hard-collinear virtualities: 

P 2 ~ QA. (50) 

Such fluctuations appear in the case of scattering collinear and soft particles and in our case they have 
momenta components which scale as 

in J- function : (pf ~ Q, pj — A, Pi± ~ y/QA), (51) 



in J'-function : (p'+ ~ A, p'r ~ Q, p' iX ~ y/QA). (52) 
Such modes are often refereed to as hard-collinear particles. 



The soft correlation function S (co,-, z/j) defined in (47) describes the contribution of the subdiagram 
with the soft momenta and low virtualities. In this particular case this is the simple product of the two 
soft propagators. Taking into account that the jet functions do not depend on the transverse momenta 
the soft part can be represented as a light-cone correlation function (CF): 

S M) K Vi )=c[ J d ±e^ J ^e— / ^ 

x (0 \q ai (Ain) q Pl (run) |0) (0| q a2 (A 2 n) qp 2 (772*1) | 0) . (53) 



In pQCD, the leading-order CF factorizes into the product of two propagators : (0|...|0) (0|...|0). 
But it is clear that this is specific for the perturbative result. In the general case one can expect general 
matrix element (0 |...| 0). Let us note that such a CF is a vacuum loop for the transverse momentum and 
simultaneously a 4-point CF function from the point of view of longitudinal subspace. Therefore integra- 
tion over transverse components is UV-divergent. Computing the convolution integrals J cLji.2 J dv\^ 
and integrating over the soft quark momenta, we reproduce the factorization breaking logarithmic con- 
tribution computed in Ref. . In our example we have only the leading-order, simple contribution from 
the hard subprocess: tree-level scattering of the transverse hard photon on the hard-collinear quark. The 
corresponding amplitude ~ 7^ can be associated with the hard coefficient function. 

The answer ( 46 1 can be interpreted in terms of a reduced diagram as in FigJ5] We observe that in 



this case the scattering process contains soft spectators and involves two large scales: the hard scale 
Q 2 and hard-collinear scale of order AQ. The presence of the soft spectators, allows us to associate the 
contribution from the soft region with the Feynman mechanism |16) . In the following, we shall refer to 
it, for simplicity, as the soft rescattering mechanism. 



B. Soft rescattering contribution for the Pauli FF Fi 

A calculation of the helicity flip FF F2 carried out in the hard rescattering picture cannot provide a 
well defined result because the convolution integral is divergent. This divergence can be understood 
as an indication that the definition of relevant regions according to Fig{T]is not complete. However, such 
a calculation allows us to define the power behavior As one can observe, F2 is suppressed as 1/Q 2 
compared to F± . This is a consequence of the helicity flip which requires us to involve one unit of orbital 
quark momentum that leads to suppression of order A/Q [one more factor Q arises from the kincmatical 



prefactor (p + p'Y in the FF definition, see Eq.(12)]. 

Consider now the diagram as shown in Fig|4J where we calculate the contribution where the hard 
photon couples to a u quark. Following the same definitions of the external momenta as before, the 
internal momenta read 

Piu=p-k 2 -k 3 , P2u = x 1 p-k 3 , p' lu =p' - k 2 - k 3 , p' 2u = yip- k 3 , 
Pig = k 2 + k 3 - xip, p 2g = h - p 3 , p 3g = yip' - k 2 - k 3 , p ig = p' 3 - k 3 . (54) 
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FIG. 4: Two-loop diagram for the helicity flip FF 



and the analytical expression for the diagram reads: 

dfadfa ^ a (p'-fa-fa)^(p- 



= C 



fah^i 



[A;| — m 2 ] 



fc. 2 — m 2 



(p 1 -fa- fa) (p~k 2 - fa) 



i'2 1 1 \ viP' - fa) T (fa + m ) i 13 (%ip - fa) 7 J 6 v-iY (fa + ™) 7 J 6 

\2 



(yip - fa) 2 (xip - fa) 2 (fa - pzY (fa - p' 3 T (fa + fa- xipY (yip' -fa- fa) 



2 ■ 



(55) 



Let us add few comments to this formula. Following conventions, we assume that the first and second 
spinor lines correspond to u quarks and we substitute instead of spinors their large components £J 2 an d 
£1^2 as defined in (30) and (31). However, we cannot perform such a substitution for all external quarks 
as we did in the case of the Dirac FF F± . In order to obtain a nontrivial helicity flip amplitude, we need 
to project the in or out collinear partonic state on the higher twist (twist-4) DAs. The projections on 
twist-4 DAs are well known and can be written in the same form as for the twist-3 case [25] . Contrary 



to the twist-3 case, the twist-4 projections do not satisfy the full set of relations (19) because twist-4 
operator includes one small component of the collinear quark field: 



twist-3 DA ~ (0| £ £ £ \p) , twist-4 DA ~ (0| £ £ r\ \p) . 
For instance, one obtains the following projector (in general, there are 9 twist-4 projections [25] ) 

t/iC <g> 75$, 



twist-4: Ty, 



(56) 



(57) 



where the u quarks projected on large components but d quark on the small component. Therefore in 
order to obtain such configuration one has to substitute instead of a d-quark spinor its small projection 
pOl: 



d(P 3 ) ^V' 3 = d(p' 3 y-f, with fj' 3 i ? 0. 



(58) 



We take into account this particular case in the expression ( 55 ) and do not consider the other configura- 



tions (with the small u-quark components) for the sake of simplicity. 

Consider first the contribution from the hard region, Eq. (23). In order to project the index \x in 



Eq. (1231) onto the longitudinal subspace, we perform a contraction 



D 



(/>) 



it m 



(59) 



Simple dimensional counting provides 
1 



1 



D(H) ~ ^ ^ Vlil ^ 2 ^3 ~ iP + P ' ] ' {n + fl} ^ ^ Tl ^ ^ 2 [Qfj * ] r3?3 ' (60) 

where we took into account the kinematical factor (p + p') ■ (n + n) ~ Q and the fact that the small 
component rj 3 is suppressed, according to Eq. ( 32 ) . Then we observe that the hard part of F 2 is suppressed 
compared to F\ (29) as expected. 
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( s) 

Consider now the soft region, expressed by Eq. (28). In the denominator D ; we obtain 



Den ~ [kl - m 2 ] [k\ - to 2 ] 

2 2 2 22 22 2 

xQ yi [(fc 2 + • «•] Q 2/i2/3 [fa ■ n] Q ii [(k 2 + fa) ■ n] Q xix% [fa ■ n] 
In the numerator, the first spinor line gives : 

t[i a {P' -fa- fa)(i + mp -fa- hh^i 

~ -2(fc 2 + fe 3 ) ' n 2/£ 7 ±£i - 2(£ 2 + fc 3 ) • n 2p to fWi £1 
Combining the contribution ~ p^7™ with the second and third lines we obtain: 

/ £> 7* (yip' ~ fa) l a (fa + m\ "f 13 (xip - fa) 7 J '6 f)' 3 f(fa + to) 7 ^ 3 

= -AmQifa ■ n) (fa + to) 7^2 V 3 P' (h + m)y&. 

The same combination of the second term ~ p' a "1^_ in Eq. (62) provides trivial results: 

P' a £ l l (W - fa) 1° (fa + to) ^ (x lP - fa) 7^2 vWCfa + m) 7 ^3 
= xi £ 2 7 l (~fa) P (fa + m)j'[p 7^2 V 3 Y(fa + m)j 3 ^ 3 
= 2xi C 2 l l (~fa) P' (fa + m) 7 { 6 fj'slHfa + to) p6 = 0. 

Therefore, we can write 

Num ~ 16 Q m(fa ■ n)(fa + fa) ■ n |» l7 J£ x (fa + to) 7^2 fj' 3 p'(fa + m) 7 ^3 



Combining Eqs. (61) and ( |65[ ) yields 
( s ) _ (p + j/) ■ (n + n) 



Q 7 



8C dfadfa 



[fc| — m 2 ] [fcf — to 2 ] 

?i7±Ci &11 (fa + m) 7 J '6 ^ n (fc 3 + m)y& 
x\x 3 1/12/3 [(fc 2 + fc 3 ) • n] [(fa + fa) ■ nf [fa ■ n) [fa ■ n] 



By simple power counting, we obtain 



w (p + p') ■ {n + n) 8 A 2 - - 

^|| ~ Q7 A ^l 1 ^2 i 2t2 %l3?3 



(p + p') • (fi + n) 



(61) 



(62) 



(63) 



(64) 
(65) 



(66) 



(67) 



One notices that we obtain the same power counting as for the hard region. Therefore we can conclude 
that the soft rescattering is also relevant for the helicity flip case and is not suppressed compared to the 
hard rescattering mechanism. 

Let us perform an interpretation of Eq. (66) in terms of hard, jet and soft functions introduced in the 
previous section. 

(s ) = (p + p') ■ (n + n) m_ 

II TO 2 



x / doJ 2 doJ 3 



hard CF 

1 [frtfL [&tfL [*&»]c 



Q 2 



yiy 3 (u 2 + w 3 )w| 



dv 2 dv 3 



1 giU 

Q 3 ^1^3 (^2 + ^ 3 ) 2 ^3 



(68) 



J'-function 



±ji.±ji (fc 2 + TO) a2fe (fcg + m) a3 p 3 Xi / + 
[A: 2 -to 2 ] [fc 2 



J-function, same as in Pi , see Eq.|(49} 

<5(w 2 - k^)8(uj 3 - fc^ )J(i/ 2 - k 2 )8(v 3 - fcg ) 



soft corr. f. S[wj,t>i] 
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We observe that the soft part and one jet function (twist-3 projection) are the same as for F 1: but the 
outgoing jet function (twist-4 projection) is different. We may expect that the soft rescattering for F2 
can also be described in terms of a reduced diagram as in Fig. [I] 



C. QCD factorization for the soft rescattering picture 



The specific feature of the soft rescattering is the presence of two subprocesses related to the two hard 
scales: a hard subprocess with typical scale of order Q 2 and a hard-collinear subprocess with typical scale 
of order AQ. Therefore description of such processes could be carried out in two steps: first, one integrates 
over hard fluctuations so that the remaining degrees of freedom describe hard-collinear and soft processes. 
From the previous analysis we may conclude that such degrees of freedom include hard-collinear ( 51|52 1 , 
collinear: 

Pc ~ (pt ~ Q, p± ~ A, P - ~ A 2 /Q) , P 2 C ~ A 2 , 

Pc ~ (p' c + ~ A 2 /Q, p x ~ A, p- ~ Q) , pf ~ A 2 , (69) 

and soft, 

p£~A, pl~A 2 . (70) 

Therefore one needs the effective theory describing the dynamics of such a system. Such effective theory, 
known as SCET, was built already for the description of heavy quark decays and some other hadronic 
reactions. Therefore we can apply it also for description of the soft rescattering mechanism. 

If Q is large enough and AQ ^> A 2 one can further use perturbation theory and factorize the hard- 
collinear fluctuations, leaving at the end only collinear and soft modes which describe soft QCD dynamics. 
Technically, such two-step factorization is described as matching of full QCD onto the soft collinear 
effective theory at the scale /i = Q (SCET/), which is equivalent to calculating the hard coefficient 
functions in front of an operator constructed from SCET/ fields described above. The second step is the 
matching of SCET/ at the scale /1 = y/ AQ to SCET//, which again corresponds to the pQCD calculation of 
hard-collinear coefficient functions (which are usually called jet functions) in front of operators constructed 
only from the collinear and soft fields. In the next section, we perform the matching of QCD to the SCET/ 
effective theory. 



III. MATCHING QCD TO SCET/ AND RESUMMATION OF LEADING LOGARITHMS 

A. Soft Collinear Effective Theory 



In this section we briefly describe the main ingredients of SCET 26-31 . The effective Lagrangian can 
be obtained from QCD Lagrangian by integrating over hard fluctuations and performing a systematical 
expansion with respect to the small dimensionless parameter A related to the large scale Q. We define 
A ~ yj A/Q, where A is the typical hadronic scale of the order of a few hundred MeV. In general, the 
physical amplitude describing a hard exclusive reaction can be defined in a convenient reference frame, 
for instance the Breit frame. Then external particles usually are hard or collinear. The fast moving 
hadron consists of energetic partons carrying collinear momentum: 

pt = (Pc • n)— +p^ c + ( Pc ■ n)— = (pt,p ±c ,p-), p\ ~ A 4 Q 2 . (71) 
The individual momentum components have the following scaling behavior: 

p£~Q(l,A 2 ,A 4 ), 



as required by (71). However, as we could see in the example above, the relevant regions could involve 
fluctuations with different momenta. We classify the different regions following the terminology suggested 
in Refs. [32 , 33 : hard ph ~ Q(l,l,l), semihard p s h ~ Q(A, A, A), hard-collinear ph c ~ Q(TA,A 2 ) or 
Phc ~ Q(A 2 , A, 1), collinear p c - Q(l, A 2 , A 4 ) or p' c - Q(A 4 , A 2 , 1) and soft p s ~ Q(A 2 , A 2 , A 2 ). 
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The large (small) components £' hc of the quark fields describing particles with momentum p' hc have 
been introduced through a decomposition of exact collinear quark fields ^>' hc '- 

e hc (x) = f < c) rf^Jx) = ^' hc , (72) 

with n £' hc — 0. The small components i]' hc are suppressed with respect to those of £' hc by a factor 
A 2 ~ A/Q and are integrated out when constructing the effective Lagrangian. 

Such definitions set the following scaling relations for the corresponding effective fields: 

e hc ~ A, n- A' hc ~ 1, A^ hc ~ A ,n- A' hc ~ A 2 , (73) 
£ - A 2 , n ■ A' c ~ 1, A c - A 2 ,n • A; ~ A 4 , (74) 

^~A 2 , g~A 3 , (75) 

with ^4'^, j4c M > ^'s^ denoting the gauge fields in the SCET, and q the soft quark field. 

After integration over hard modes we reduce full QCD to the SCET/ which describes the interaction 
of particles with hard-collinear and soft momenta. This theory still includes the particles with large 
virtuality of order AQ 3> A 2 if Q is large enough. Therefore, if possible, one can perform a matching of 
the SCET/ to the effective theory which contains only collinear and soft particles (SCET//). In present 
paper we consider in detail the matching of QCD to SCET/ and resummation of Sudakov logarithms 
which arises due to the evolution of the SCET operators. 

The effective action describing the interaction of the hard-collinear and soft particles can be written as 
an expansion with respect to A [30], [31] 2 : 



where 



4 0) = J d 4 x & c (a;) [in-D + gn- A s (x.) + W x ^ it> ± ^ ^ hc {x), (77) 

4? = j d ' x + q{x-)W'\x)iD^' hc {x)- (78) 

where X- = ±{x ■ n)n. The fields {£' hc , A' hc } and {q 7 A s } describe hard-collinear {(p' hc ■ n) ~ Q) and soft 
fields respectively, whereas the covariant derivative reads iD = id + gA' hc . The Wilson lines are defined 
as 



W'(x) = Pexp i^ig j ds n ■ A' hc {x + sfi) j . (79) 
We do not write explicitly the gluon part £ym because we will not use it in this work. C s is the usual 



QCD Lagrangian with the soft fields. In Eqs. (77) and |78|) we provide only the contributions which will 
be relevant for our discussion. 



The expressions in Eqs. (77) and (78) have definite (homogeneous) scaling in A which is indicated by 
the number in the superscript brackets. In order to achieve this the arguments of th e so ft fields are 
expanded with respect to the small components of the position arguments 3 . From Eq. (77) one can see 



2 There are two different technical formulations of SCET developed in 27 29 -31 . In the present paper we follow the 
technique suggested Bcneke et al. in Ref. [3D]. 

3 In momentum space such "multipole " expansion corresponds to the expansion with respect to small momentum compo- 
nents. 
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that the soft-gluon fields couple to the hard collinear fields only via the longitudinal component n ■ A s . 
Using the field redefinition 

i' hc {x) -> S n (x.)^(x), A' hc {x) -+ S n {x-)4°\x)Sl(x-), (80) 
with the soft Wilson line 



S n (x) = P exp sig ds n ■ A s 

I J — oo 



(x + sn) \ , (81) 



we can eliminate the soft field from the leading-order Lagrangian (77). However the soft Wilson lines S n 
remain in the external operators with soft fields in order to ensure the gauge invariance. 

Obviously, all the above results also hold for the second collinear region with momentum (phc ■ n) ~ Q, 
by merely interchanging light-cone vectors n O n and substituting corresponding hard collinear fields 
{6io A hc }. 

The formulation of SCET described above can be extended by introducing the so-called soft-collinear 
or messenger modes as discussed in [33]. However, such particles have virtualities which are much smaller 
then the typical hadronic scale p 2 sc <C A 2 . This situation was investigated in detail in several papers ( 
see e.g. [35H37]). It was shown that the existence of such modes depends on the precise form of the IR 
regularization used in massless pQCD. Therefore it was suggested that in the processes with real hadrons, 
where all nonperturbative effects have typical scales of order A, such low-mass degrees of freedom cannot 
appear because they are clearly an artefact of perturbation theory. Therefore we do not include them in 
the present considerations. 



B. Construction of the operator basis and leading order coefficient functions 

In this section we briefly describe the matching of QCD to the relevant leading-order operators in 
the SCET/. The leading-order matching the e.m. current onto the SCET operators has been already 
introduced and studied earlier. In Ref. [35] it was used for a description of DIS at large x — > 1 and in 
Ref. [39] for the description of Drell-Yan production. The matching onto subleading operators was also 
discussed in [45] . For the convenience of the reader we repeat here the main steps of these calculations 
in order to introduce required notations. 

In order to obtain the allowed SCET operators we take into account the restrictions imposed by the 
SCET counting rules, gauge invariance, and invariance under the reparametrization transformations. 
Explicit construction of such operators can be performed in the same way as it was done for heavy- 
to-light transitions in the works of Refs. [4"0T|44| . The building blocks, invariant under collinear gauge 
transformations are well known and read 

{(ZLW), (W^ hc )}~\, (82) 
{[W'UD^W], [Wt*D M W]} = {^,^}~A. (83) 

In the terms with [...], the derivative is only applied inside the brackets. 

For the LO operator one can easily construct the expression which consists of two quark jets: 

0(8 lt 82) = [e hc W] (sin) ® [wHho] (s 2 n) (84) 

where we do not explicitly write the color and spinor indices for each jet and the symbol ® is used to 
stress that their indices are not contracted. From the previous discussion it is clear that such an operator 
is relevant for the Dirac FF F% . For the case of Pauli FF F% we need the subleading operator involving 



the transverse gluon field as in Eq. ( 83 ) . In this case it is useful to take into account the constraints 
imposed by reparametrization transformations. The details were already discussed in the literature, and 
we refer to Refs. [30] HQ1 E3 E3] for these details. The relevant for our consideration subleading operators 
can be written as 

On(si, s 2 , a 3 ) = WhcW] (s 2 n) / x (s 3 n) ® [W^ hc ] [am), (85) 



) = [i'hcW] (s 2 n) ® 4±M [W^hc] (sin). (86) 
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So that for matching of vector current we can write 



g(0) 7 ^(0) = J dhdh tr [C»{s x ,h) 0«( Sl ,s 2 ) 

+ I dhdhdh tr C£{h,h,h) Ol(sx,s 2 ,s s ) + C{ l l (s 1 ,s 2 ,s 3 ) 0£(si,s 2 ,s 3 ) (87) 



where Sj = s^Q. The coefficient functions C arc defined as matrices in the spinor and color indices and 
the trace has to be understood in a sense of contractions of all the indices between coefficient functions 
and operators, for instance: 

tr [C»{a x ,h,Q/n) 0( Sl ,s 2 )} = f C»(h, h)] [& c a W] ( Sl n) [W^ hc p ] (s 2 n). (88) 

J aft 

The further details of our calculations arc presented in Appendix A. It is convenient to pass in mo- 
mentum space where the final result can be presented in the compact form: 



g(0) 7 M 9(0) = C A (Q, M ) 



dr C b (t,Q,h) O b [t] + 



(89) 



where the scalar coefficient functions Ca.b include all relevant contributions with large logarithms, and 
the operators are defined as 



07 = (?HcW) (0) 7 ^ (W^ hc ) (0) = (eW) 7 ^ , 



B [r] = ^ he W') (0) 



2tt 



-is(P -n) r Af 



4' ± (sn) + e"( p -") (W+ffcc) (0) 



(90) 



(91) 



From the tree-level calculations it follows that 

C A (Q^ = Q) = l + 0{a s ), C B (T,Q,tx = Q) = l + 0(a s ) . (92) 
Note that the SCET operators depend also on the renormalization scale fj, that was ignored for simplicity. 



C. Resummation of large logarithms 

The next important step is the resummation of the large logarithms or, equivalently, the solution for 
the evolution of the SCET/ operators. As we described above, we expect that the scale for the remaining 
hard-collinear subprocesses is of order AQ. Therefore it is natural to set the factorization scale /j, 2 to be 
of order AQ. However, we then obtain in pQCD large logarithms InQ 2 //! 2 which must be resummed to 
all orders. Such resummation can be easily performed with the help of the renormalization group (RG) 
and has been carried out for many applications. We therefore only briefly describe the main steps and 
provide the final results. 



We start our discussion from the coefficient function Ca in front of the LO operator, Eq. (90). The 
corresponding RG equation reads 



dlnfj, 



C A (Q, M) = C A (Q, h) 7 (a.) , C A (Q, h = Q) = 1 + 0{a s ), 



The anomalous dimension 7 (a s ) is defined by renormalization of the operator Oa ( see, e.g. 
well known that to all orders the anomalous dimension 7 can be represented as 



(93) 
]). It is 



7(««) 



(a s )\n^- +7(a s ), 



(94) 



where the coefficient in front of the logarithm in Eq. ( 94 1 is known as the universal cusp anomalous 
dimension, and controls the leading Sudakov double logarithms. Such specific term is usual when Sudakov 
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logarithms appear for the quantity under consideration. The single-logarithmic evolution is controlled by 
the j(a s ). 

The solution of Eq. ( [93] ) provides a systematic resummation of large logarithms in pQCD. In order to 
find Ca in the next-to-leading logarithmic (NLL) approximation, 



C^ LL ~ exp \J2 a ™< ln " +1 + & « a " ln "} » 
one needs to know the 2-loop cusp anomalous dimension r .i and the leading-order term 7^ 



Tcusp (^s) 



47T 



47T 



where 



r = 4C F , Ti = AC f 



67 tt 2 \ 10 



) 71 



-6C, 



F ■ 



The explicit NLL solution reads 

CT L (Q^) = e- s ^"^ U A (ji h ,n) [l + 0(aM) 

where 



S{Q,Hh,n) = -^lnrln^ + 

Po v ^Po 



47T 



In r — 1 



1 



Pi , 2 



ri _ ft 
r A) 



— 1 — In / 



(95) 



(96) 



(97) 



(98) 



(99) 



U A (fi h ,fi) =r 2/3 o . 
with r = a s (n)/a s (/j,h) > 1 and /?— function coefficients 
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: C A --n f , ft 



34 „ 2 

t 6a 
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(100) 



(101) 



In Eq. ( 98 ) we assume that evolution is running from the initial scale fXh (which should be of order Q) to 
scale fj, of order (QA) 1 / 2 . 



A similar technique can also be used for the subleading operator of Eq. ( |9l| ) . Notice that in this case 
our calculation also provides the practical check for the existence of the convolution integral in Eq. ( p9| . 
If it does not exist then our suggestion about the factorization must be reconsidered. 

In order to find the anomalous dimension one has to compute the diagrams shown in Fig(5] The 




FIG. 5: One-loop diagrams required for renormalization of the three-particle SCET/ operator. As before, soft 
gluon lines indicated by an index s. The wave function renormalization diagram is not shown for simplicity. 



operator RG equation reads 



- J cIt'j[t,t'} O r [t'} 
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with the evolution kernel j[t, t']: 

i[t,t'] = j- 5(t - r')r cusp [a s ] In 

where (c.f. [15]) 



Q 2 



y[r,r']= -S(r-T r ) (c F 


"5 . _" 
--lnr 


c A f\ c A ■ 


+ (<v-£) [ 


% 0(r'<f)+0(f <t')- 

T T 


+ 



(r < t') <9(t' < t) 



(102) 



(V - r) (r - t') 
C F f 



2 



3 

f ~ 2 



+ fl( T ' <T ) 

T \ Z T r 



(103) 



where the prescription [...] + is defined for the symmetrical kernel /(t,t') = /(r',r) as 

[/(r,r')] + = | r'/(r,r') - 0(r)] . 

Computing the convolution integral with the LO Cb[t] = 1 yields the well defined expression: 



/ dT , V[ T ',T} = - 
Jo 



2C p — t;Ca 



(104) 



which does not depend on r. Hence we can conclude that the leading logarithmic convolution integral in 
Eq. ( 89 1 is also well defined. 



The corresponding RG equation for the coefficient function reads: 



H-^C B (T,Q;fj) = -r CU8 p[a a ]ln ( j-: j C'i;{ t. (): // i 



^ f dr 1 V[t',t}C b (t',Q;m) 

T JO 



(105) 



Similar equations have been studied already in heavy- light decays (see, e.g., Refs. 
solution of this equation can be written as 

C£ ll (t, Q; m ) = e"^'") / rfr'C/ [r, r'; ^, /*] C^ 0) (t',Q; 



where the evolution kernel satisfies the integro-differential equation 

n 1,1 
r/ Ft. t : ii.u. n] — — / 



M ^[r,r';^,M] = — / rfr" V[r", r]J7 [r", r'; /i/j, /z] 



The NLL 
(106) 

(107) 



with initial condition {/ [r, t'\ \ih, fJ>h] = 5(t — t'). Recall that, in order to sum the large logarithms 
the initial scale \Xh should be of order Q, and the evolution ends at /i of order (AQ) 1 / 2 . The Sudakov 
factor S(Q; fih:^) 1S the same as in Eq. ( |99[ ). Taking into account that at NLL approximation the initial 
condition Cr (t', Q; fJ,h) is given by the tree-level expression Eq. (92 1, one can perform the integration 
over r' in Eq.(106), yielding : 



with 



d a f 1 

H-tUb [r;/x fc ,/x] = — / dr' V[t',t]U b [t'; /j, h , 

dfi 7T J Q 



(108) 
(109) 



and U b [t; fih, A*/i] = 1- The solution of this equation can be found numerically. We have found that to a 
very good accuracy the approximate solution can be written as 



U B [t; fi h ,fi] « \jj, h , fj] = 



/ ln[^/A(»/) ]^ 27B//3 ° 
V ln[^/A("/)] 



(110) 
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FIG. 6: Left panel: numerical evaluation of the evolution kernel Ub[t,Q, y/QA] as a function of r at fixed Q 2 . 
Right panel: Ub[t, Q, \/QA] as a function of Q at fixed r = 0.5 for two different values of A (solid curves). The 
approximate solution U B PP [/J.h,fj] of Eq. (110 1 is shown by dotted (blue) curves. 



used for calculating 



with effective anomalous dimension 7# defined in Eq.(104) and with soft scale A^ r 
the running coupling a s . To avoid confusion let us note that the soft scale A which we used to define the 
hard-collinear scale ~ QA is diff erent A ^ A^ nf \ This difference provides the slow dependence on Q in 
the approximate solution of Eq. (110). In Fig. [6] we show Ub [t ; Q, VQA] computed for different values 
of Q and A and compare it with the approximate solution [/^ pp . We obtained that for all considered 
cases to a very good accuracy the kernel does not depend on the momentum fraction r and evolves quite 
slowly with respect to Q according to (110). At the end let us note that the similar approximate solution 

Ha. 



for the single-logarithmic evolution was also found for the heavy-light current in Ref. 

The obtained results already lead to some qualitative features when applying this formalism to the 
proton FFs, as will be discussed in the next section. 



IV. QCD FACTORIZATION AT LEADING ORDER USING THE SCET APPROACH 



In this section, we consider the matching on SCET// and discuss the factorization formula for the 
soft rescattering mechanism. We perform an analysis of the dominant regions using the methods of the 
effective theory. We restrict our consideration only to the terms relevant at leading logarithmic approxi- 
mation both at SCET/ and SCET// levels. The general, all order analysis is much more complicated and 
goes beyond our present considerations. However, using the results obtained above, we suggest a leading 
order factorization formula (i.e. restricted by leading logarithms) for the Dirac FF Fi which includes soft 
and hard rescattering contributions. 

In this section we would like to demonstrate that the soft rescattering contribution can be estimated 
in SCET using the counting rules ( 73||75 1 without direct calculation of the diagrams as we did before. 
Such counting is an important ingredient of a factorization proof and can be considered as a quite general 
argument in support of the nontrivial soft rescattering contribution. 

In addition to the field relations, we also need the counting of the energetic (collinear) hadronic state. 
It reads 



\Pc) ~ A" 



(111) 



and follows from the conventional normalization ( 13 1 



Let us start from the well known hard rescattering picture. From existing results one can easily obtain: 
(p' \J»\p) N'~/1_N F t ~ N'l^N 1 -^ ~ N'^N A 8 , (112) 

and 



,"■ A ) = {P + PT N'lN F 2 (Q 2 ) 



m N 



m N Q 



(113) 
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where the nonperturbative scale is presented by the nucleoli mass m^r ~ A and the overall normalization 
constant for the nucleon distribution amplitudc//v ~ A 2 . 

The same counting in SCET is obtained directly from the dimensional analysis of the leading op- 
erators constructed from the collinear and soft fields, which represent the main degrees of freedom of 
SCET//. However, difficulties arise due to nonlocal contributions with inverse powers of 1/A momenta 
if one computes time-ordered products involving SCET/ fields. We shall follow the strategy suggested in 
Refs. [mill]. Then matching onto SCET//, 

(114) 



SCET/ [p 2 hc ~ QA, k 2 s ~ A 2 ] "^X SCET// [p 2 ~ A 2 , k 2 s ~ A 2 ] 



can be performed in two steps: decoupling the soft fields from the hard collinear modes using field 
redefinitions expressed by Eq. (80), and subsequently matching hard-collinear modes to collinear ones 



rc'(0) a'W a") a(")\ , id a' p 4 1 

\Shc ' A hc >?ftc' he I i^c' A c^c, sich 



,'(o) AO) A (o) 



(115) 



lowering the off-shellness of the external hard-collinear fields. Notice that the last step changes the power 
counting of the fields from Eq. ( 73 ) to Eq. ( 74 ) . 



In the case of the hard rescattering we perform the matching of QCD directly onto SCET//. Therefore, 
the power counting is simple because it does not involve the intermediate effective theory. 

Matching for the Dirac FF F\ involves the six-quark operator constructed only from the collinear fields 
£' c , £ c and Wilson lines with longitudinal collinear gluons n ■ A' and n ■ A respectively. It is the product 
of two twist-3 3-quark operators which define the leading twist nucleon DA (14). Then using ( [74] ) one 
obtains 



{p' \j£\p) ih) ~ (p' \&W0)*Cl(Q)*(0\ZeUc\p) ~ A 8 NY^N. 



(116) 



For the helicity flip FF F2, the matching involves the product of twist-3 and twist-4 operators, as we 
discussed in Sec. Ill Schematically this situation can be described by substituting £ c — > rj c ~ £ C /Q. Then 



P 





\(h) 




TP 


p) - 





p' e c &' c \0)*q(Q)*(0\r, c Uc P) - (n 



-nY A 10 N'lN 



(117) 
(118) 



and we obtain that F2 is suppressed as 1/Q 2 relative to Fx as it should be. 

In order to estimate the soft rescattering contribution one has to perform a more complicated analysis 
with the two-step matching: QCD— >SCET/ — >SCET//. The matching onto SCET/ has been done 
in Section IV and for the electromagnetic current at leading order it yields the formula Eq. ( 89 1 . In 



matching onto SCET// we need at least six collinear quarks in order to have overlap with the in and out 
nucleon states. Next, guided by the perturbative QCD calculations from Sec. [n]we need higher order 
vertices in order to describe the soft spectators in the intermediate state. 



A. Leading-order SCET analysis for Dirac FF Fi 

Let us start with a discussion for the FF F\ . To leading order in 1 /Q we can restrict our consideration 
by the first term in Eq. ( 89 ) . Therefore our task is to compute the time-ordered product which contributes 
to the matrix element 



(p'\Jl (0)\p) {s) ^( P '\ T{C a O 



(n) 

SCET7 



-tC, 



\p) 



(119) 



The calculations amount to integrating out hard-collinear modes and, if possible, to deriving the expres- 
sion for the vector current J±_ in terms of SCET// collinear and soft fields which can be schematically 
written as 



T { C A Ofr • 



SCET, 



+iC. 



C A Tr 



IT 



'{Oout(v'c) e l£ ^}*J'*T{S(^) e i£ *}*J*r{0 in (O e^"'} 



(120) 
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where J and J' are jet functions, C c n ^ denotes the collinear Lagrangian, Tr denotes contractions over the 
Dirac and color indices which are not shown explicitly, and where we used the notation (p c>s = {£ c , s , A cs }. 
We also assumed that the collinear operators Oi n ,out have nontrivial overlap with nucleon states. The 
Dirac matrix 7^ is associate d wit h the vertex of the Oa operator. It is clear that the existence of 
the factorized representation (120) is equivalent to establishing the factorization theorem. Guided by 



our QCD analysis, carried out in Sec. II, we demonstrate below that at leading order in 1/Q such a 
contribution definitely exists. For simplicity, we restrict our consideration to a leading-order analysis in 
a s , and consider it as a first step towards a complete proof. 



Obviously, the time-ordered product in left hand side of ( 120 ) can be represented as the product of 
two: 



T{...}=T (? hc W(Q)e 



iC 



SCET 1 



71 T W£ hc (0) 



iC 



(n) 

SCET, 



Tout ^fj_ Tit 



(121) 



where we "freeze" the soft fields, i.e., consider them as external. As calculations of the each of the T 
products are almost identical, we only consider one of them. The result of the integration over hard- 
collinear modes can be schematically written as 



Tout = T £ hc W e^soBx, 



(122) 



where the last equation shows the desired structure in terms of collinear and soft fields. Combining such 



results for T in and T out we obtain desired representation (120). 



Let us now consider in details the calculation of the right-hand side of Eq.( 122 ). The relevant T product 
is of order A 3 and to leading order in a s reads 



= T 



(3) 



d A X\ 



d i X2 



d 4 x 3 



T 



-(i) 

-«3 



(123) 



where C^ q is the leading-order soft-collinear contribution in Eq. (78). We did not find the other possibil- 



ities to obtain the leading in the 1/Q result. Time-ordered products with insertions of other higher order 

(n) 

contributions C^ q with n > 2 from the collinear or soft-collinear sectors can provide only suppressed 
operators in SCET// and therefore can be excluded from the consideration. Performing a decoupling of 
the soft field with the help of Eq. (l80|) we obtain 



eL^'(o)^e;i V(Q)st(o) 



(124) 



The eikonal factors ensure the gauge invariance of the soft sector described by the soft quark fields. 

Subsequently, we compute the contractions of the hard-collinear fields which can be conveniently presented 

(3) 

by Feynman graphs. The leading-order contribution to T^J t is given by the set of diagrams shown in 
FigjTJ Note that the last two diagrams with the three-gluon vertex have zero color factor and therefore do 





FIG. 7: Leading-order SCET diagrams required for the calculation of jet functions. The inner dashed and curly 
lines denote hard-collinear quarks and gluons, external dashed lines correspond to collinear quarks, fermion lines 
with crosses denote soft quarks. Black squares denote the vertex of the SCET/ operator. 

not contribute. This is in full agreement with the similar observation made in Ref. [T^]. The remaining 
diagrams have a similar topology and the corresponding power counting can be easily established. The 
contractions of the hard-collinear fields yield 

((ft • A hc ){x 2 )(n ■ A hc )(x 3 )) - A~ 2 , (125) 
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(126) 



i.e., all hard-collinear contractions cost A -2 , which results from the hard-collinear propagators in momen- 
tum space. Remember that we assume that external hard-collinear particles are matched onto collinear 
ones. Therefore taking into account the external collinear and soft fields we obtain 



T. 



(3) 



3 \3 



A A 



x A 2 A 2 A 2 



A 



(127) 



h-coll contractions 2 soft fields 3 coll fields 



The same counting is also relevant for the second time-ordered product in Eq. ( 121 ). Therefore the order 
of the total contribution in SCET77 now reads 



(P'\ J1\P) 



T (3) y T (3) A -4 
1 out X 1 in X A 



N'^N ~ A 8 N'jIN. 



(128) 



We observe that the soft rescattering contribution has the same power suppression as the hard one [ 
Eq. pl6| ]. 

Let us briefly discuss the general structure of the soft rescattering contribution. It is clear from the 
above consideration that the leading-order jet functions can be computed from the diagrams in Fig. [7] 
The details of their calculations and explicit expressions will be presented in a different publication [50f . 
From the QCD calculation, we noticed that at tree level the transverse momentum is completely defined 
by external soft and collinear fields and therefore it scales as k± ~ QA 2 . Such counting ratio remains true 
for the hard-collinear lines inside diagrams due to the momentum conservation. Therefore the transverse 
components in the hard-collinear propagators (for tree diagrams only!) can be neglected. Consequently, 
the arguments of the external collinear and soft fields are local in transverse space 4 and depend only on 
the relevant light-cone components. In SCET the same properties follow from the multipole expansion of 
the fields with respect to the small parameter A. Therefore computing the diagrams in Fig. [7] and passing 
to momentum space one obtains 



T out- J Dyi0 out (yi) J rfwi :2 J' (y t ,uj t Q) S n (ui), 



(129) 



with the following collinear and soft operators: 

3 



out {y i )=A Y[ 



dz? 
2ir 



e ~\ { p '- n )(yi z i +y2*i +V3*t) 

xe^ k OY?(izf n) ? c W?(h4 n) ^W' c \h4 n), 



(130) 



and 



J'j'k' 



dz 1 



2tt 



^gj( w l 2 l +"2^ 2 ) 



k'l 



k' 



(131) 



where we do not show for simplicity the spinor indices, P' denotes the total collinear momentum operator, 
= {zi ■ n), zf = (zj • n), and dz^~ = dz+Q. The structure for can be obtained in an analogous 



way. Combining these results we obtain operator expression with structure ( 120 1 which schematically 
can be written as 



T{C a O^} ~ C A Tr [^0 ou t(yi) * J' Ui Q] * {S R S„} * J [a*, u t Q] * O in (a 



(132) 



Substituting these results into the matrix element of Eq. (119) and taking account of the decoupling of 
the collinear and soft modes (see, e.g., [23111]) one obtains three matrix elements: two with the collinear 



4 We choose x = in the ( 119 I that correspond to xj_ = 0. 
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fields and the soft correlation function. The collincar matrix elements can be easily converted into DAs 



DVi O out 0i) |0) J'(y i ,w i Q)= / Dy, J'G/^Q) 



(133) 

Rewriting the initial matrix element (0| £ c £c£c \p) ~ ^ in the same way and combining all contributions, 
we obtain the factorization formula for the soft rescattering contribution: 

x J Dxi^(xi,fiu) J dv x dv 2 3(xi,ViQ,fj,i,fj,n)S(iJi,Ui',fj,ii), 



(134) 



where the soft correlation function is defined as 



a( \ f dm f dr] 2 

with the operator 



' / ^e lAlWl+jA2 " 2 (0|O s (^,A,)|0), (135) 

27T / 27T 



1 J 



fc' 



13 V 



(136) 



which is shown graphically in Fig. [8} In the last equation we assume qSn(ji2n) = q{r]2n)Sn(r]2n), and the 

r 



f 
k' 



ii ^ i 

.i' i j 

ii' -fc 



FIG. 8: Graphical representation of the operator of Eq. ( 136 1 entering the soft correlation function. Dashed and 
dotted lines denote the different light-cone directions. 



color and Dirac indices are shown explicitly. Furthermore in Eq. (134), Ca denotes the hard coefficient 



function which has been computed in the leading-order approximation ( 92 ) . 

In Eq. (134) we show explicitly two factorization scales fii and fin . The total contribution, as usually, 



does not depend on these auxiliary quantities. The first scale \ii arises at the matching QCD to SCET/. 
The evolution equations at leading logarithmic approximation with respect to \xi were discussed above. 
In practical applications it is convenient to fix this scale at the value ~ AQ. Then the large logarithms 
ln[<5 2 //ij] can be resummed solving RG equations. The second scale fin appears when one performs the 
second reduction to SCET77. Usually, the corresponding coefficient functions ( = jet functions J' and J ) 
are computed at fi 2 u ~ AQ and then the scale is fixed to be fi H ~ A 2 . Again, arising large logarithms 
ln[QA//i|j] must be resummed with the help of evolution equations for nucleon DAs ^{x^fin) and CF 
S(uJi,Vi, fin). The evolution of the DAs is well studied in the literature (see, e.g. Refs. [HUE] and [ST] 
for recent progress ) but the corresponding equation for S(uii,Ui, fin) is new and has not been derived 
before. Such a calculation must be done because it will provide an important check of the factorization 
formula (134) at leading logarithmic accuracy. A derivation of the jet functions and evolution kernel 
for S(uii,Vi, fi) will be presented in a separate publication |50j . In the Appendix B to this paper we 
demonstrate how the perturbative QCD result of Eq. (48) is reproduced from the corresponding SCET 
diagram in FigjT] 
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It turns out that the product of the nucleon DA and jet function has the same Dirac and color structure 
as pi: 



+ (l A (xi,Wi) P+[^ 5 C} ap [N+] a 
+ n T {x h Ui) P+ [tfaxClnp [7^75^+] 



(137) 



where the coefficients Qx{Xi,0Ji) are linear combinations of the nucleon DAs (16) and hard-collinear jet 
functions 



tt x (Xi,Ui) = Jxv{xi,uJi)V(xi) + J X A{%i,Ui)A(xi) + Jxr(x l ,uJi)T(x l ). 



(138) 



Therefore, the jet function can be interpreted as a hard-collinear component of the three-quark nucleon 
wave function describing the transition of the three collinear quark state into configuration with one hard- 
collinear and two soft quarks. Correspondingly, the CF S(u3i,Vi, /zj/) (135) describes the propagation of 
the soft diquark state in the background of the soft-gluon field created by a fast moving active quark: i.e., 
it describes the soft overlap of the nucleon wave function. Therefore we expect that the soft rescattering 
picture can also be associated with the well known mechanism, suggested by Feynman a long time ago 



In the factorization formula of Eq. ( 134 ) we restricted the fractions u>i and Vi to be defined on the real 



semiaxis, assuming that all the functions in (134) are real functions. This allows us to avoid the poles in 
the propagators of the tree diagrams in Fig. [7j and ensures that the jet functions are real. Recall, that 
the reality of the nucleon form factors is guaranteed by the time reversal invariance of QCD. 

The other interesting observation which follows already from the QCD computation ( 48 ) is the absence 
of the end-points singularities in the convolution integrals of DAs with jet functions in Eq. (137). Let 
us assume that the convolution integrals with respect to the soft fractions u>i,Vi in (134) are also well 
defined. Then this allows us to suggest that hard rescattering and the soft rescattering mechanisms 
provide additive contributions to the total FF Fi at least to leading logarithmic accuracy: 



F\ (s) + F} 



CO 



(139) 



with the well known expression for the hard rescattering part: F^ = \& * H * 

convolution integrals for F± at leading logarithmic accuracy are also well defined 
conclu de th at there is no double counting in this case. The formula of Eq. ( 139 1, together with the result 
of Eq. ( |134[ ), is our suggestion for the full factorization formula for the Dirac FF at large Q 2 . We would 
like to emphasize that the obtained results have been derived at leading order and only partially verified 
at the leading logarithmic approximation. A discussion of an all order factorization proof for Eqs. (134) 
and ( 139 1 requires more detailed analysis and goes beyond this publication. 



Recall that the 
Hence we may 



B. Leading order SCET analysis for Pauli FF F 2 

In contrast to Fi, the description of the Pauli FF F^ is more complicated. First, it is well known that 
hard gluon exchange can produce large logarithms, 

F^ h) = **H**~a 2 ln 2 Q 2 /M 2 , (140) 

which arise due to the end-point singularities in the convolution integrals |52| . This is perhaps, an 
indication that the hard and soft rescattering mechanisms overlap. Therefore, in order to find the correct 
description for F% one has to formulate a clear recipe for how to avoid double counting in the calculation 
of soft and hard rescattering contributions. Such a problem for F 2 , probably, arises already at the level of 
matching QCD to SCET/. However, our analysis from the previous section does not show any problems 



with the SCET/ convolution integrals for the coefficient function Cb [see, e.g., Eq. (104)]. Moreover, 



resummation of the leading Sudakov logarithms can be carried out exactly because the problematic 



logarithms ( 140 1 admix only at the next-to-leading accuracy. The structure of the logarithms beyond the 



leading order is an important subject which remains to be established for a full proof of the factorization 



24 



theorem. As a first step in this direction, we demonstrate here that the SCET counting rules confirm 
the power suppression of the soft rescattering contribution in F 2 obtained from the QCD calculation in 
Sec. M 



For this purpose, the relevant part of the SCET/ vector current follows from Eq. (89) as 



4 (o) 



i 
Q 



(141) 



with ^.j_( T ) + -A-l( t ) as given in Eq. (91 1 . For the qualitative discussion we consider the first term with 
only. Following the similar arguments as before we arrive at an analysis of time-ordered products: 



rp{i) _ T 

-'-out 



T / T»((\\ J^ n) +iC^\ _ 71 T (4) T (3) 

1 \ J \\ [V > e J - (fi-d) out 171 

{M (0)A(0)e l£<S) } , r£> =T{(WHhc) (0)e i£W } 



(142) 
(143) 



term 



(3) 

The second term Tr' in Eq.(143l is the same as the one appearing in F\ Hence we only 



To obtain Eq.(142) we converted (141) into position space and for simplicity wrote explicitly only the 



have to consider the new term T^ t . Consider the following contribution: 



^l^lf:(0)4i(0) / d Xl C^{ Xl ) I dx 2 L^{x 2 ) I dx 3 C^(x 3 ) 



'(0)/ 



(144) 



~ T (t> c (0)Acx(Q)Jd Xl ? hc A hcX q(x 1 ),fdx 2 g (n ■ A hc )£ hc (x 2 ) J dx 3 ff c (n ■ A hc ) q(x 3 j) , (145) 
where we substituted the small component 

ff c = £^± c (m • ^c)- 1 ! ~ £ gA c± (in ■ ty't (146) 



The presence of the small component fj' c can be explained by interaction with the longitudinal photon. 
The outgoing collinear state must have one collinear transverse gluon or transverse derivative in order 
to satisfy conservation of the orbital momentum. Contracting the hard collinear fields in Eq. ( 145 ) one 
obtains the diagram as in Fig. [9] Using Eqs. ( 125 ) and ( 126 ) one easily obtains 




FIG. 9: SCET diagram for matching the quark-gluon hard-collinear state onto three collinear quarks for Pauli 
FF. The small component fj' c is shown as a quark gluon state £' C A C ±. 



T, 



(4) 



A 3 A 3 



2 soft filcds 



aW 

3 coll fields 



A~ 2 A _2 A~ 2 



h-coll contractions 



Then the total contribution reads 



To ( «l x if > x A- 4 (n + nTNiN 

A 6 A 8 A~ 4 (n + n)» N IN ~ A 10 (n + N IN, 



(147) 

(148) 
(149) 



i.e., we obtained the same result as in the case of the hard rescattering mechanism ( 118 ). In the Appendix 
we demonstrate that the diagram in Fig. [^reproduces the QCD expression for J' from Eq. (68). However, 
contrary to the Dirac FF, the convolution integral with respect to the collinear momentum fraction * J' 
is not defined due to the end-point divergencies. Therefore we assume that the matching onto SCET// 
for the Pauli FF F 2 cannot provide a well defined expression. As we discussed above, even matching onto 
SCET/, one is faced with the mixing problem between hard and soft rescattering contributions. 
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V. PHENOMENOLOGICAL APPLICATION TO THE NUCLEON FFS 



In order to perform a first phenomenological analysis we introduce SCET/ form factors denned as the 
following nucleon matrix elements, 



(p'\(? q W')7l{wH q ) \p) 



SCETj 



N(p') fll fN(p) /f(Q J)M ) = N' + y»N + f?(Qj^), (150) 



and 



(p\Ob[t] \p) SCETi = N' + 1 N+(p) !^*(r;Q />A i) 



(151) 



with the operator Ob[t] dehned in Eq. (91 1. We also used the notation Qi ~ y/QA in order to stress that 
the dehned quantities do not depend on the large scale Q 2 . We indicate explicitly in the right-hand side 



of Eqs. ( 150|15l' l the renormalization scale dependence. Taking the nucleon matrix element from both 





FIG. 10: Graphical representation of the SCET/ FFs /i (left panel) and /2 (right panel). In /2, r corresponds 
with the gluon momentum fraction. 



sides of Eq. (89 1, we obtain 

F 1 (s) (Q 2 ) = C A (Q, W )^e 9 f!(Q I , t x) = C A (Q,fi) /i(Q/,/i), 



(152) 



F 2 (s) (Q 2 ) 



Q 2 



" l N 

Q 2 



dTC B {T,Q,») / 2 (r;Qj, M ). (153) 



Eqs.(152| and (153) are presented in Fig 10 in graphical form 



Using NLL approximation for the coefficient functions (98) and (108) these results can be represented 



as 



F^(Q 2 ) 



t iV ( ,-S(Q;Mh,M/) 



Q 2 



/ dTU B [T -^h,^} f 2 (T- 1 Q I: fl I ) 

Jo 



(154) 
(155) 



where the scale [ii ~ Qi ~ y/AQ. From the right-hand side of Eqs. ( 154|l55 1 one can see that the SCETi 
FFs /i t 2 depend now only on the hard-collinear scales. All dependence from the large scale of order Q 2 
is factorized into Sudakov factors e~ s ^' fJ,h,,1I \ This is the main feature of the Feynman mechanism. We 
could expect that the hard scattering contribution in F\ provides corrections of order a™ +2 In™ Q 2 //i 2 
which are suppressed relative to the contributions computed in Eq. ( 154 ), and therefore can be neglected: 



Fl {Q 2 )^F[ s \Q 2 



(156) 



In the case of the Pauli FF F 2 , the situation is more delicate due to a possible overlap of hard and soft 
rescattering terms. From the calculations of the hard scattering contribution [S5], one obtains, due to 
the end-point singularities, contribution of order a 2 In 2 Q 2 / /i 2 . Such logarithms are of the same accuracy 
as next-to-leading Sudakov or single logarithms in Eq. (155). Therefore, Eq. (155) is exact only at the 



level of leading Sudakov logarithms. Beyond this accuracy one has to perform a more accurate analysis 
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in order to avoid double counting. For a first numerical estimate, we shall neglect the hard scattering 
contribution in F 2 assuming 



F 2 (Q 2 ) 



F 2 (s) (Q 2 



(157) 



Such an approximation, perhaps, may work if Q 2 is not very large, of order a few GeV 2 , and one may 
expect that the dominant contribution is provided by the soft spectator contributions of Eqs. (154) and 
( fl55| ). 

First, it is inter estin g to investigate how strong suppression is obtained from the resummed Sudakov 



logarithms. In Fig 11 we demonstrate the results for the leading Sudakov logarithm factor e s (Q^^,^i) 




FIG. 11: Effect on the FFs from the leading logarithmic summation (LL) for different values of A. 



taking /i/j = Q and fij = ^AQ. We use two different values for the soft scale A = {0.4, 0.6}GeV and 
consider Q 2 = 4 — 16 GeV 2 . For our numerical estimate, we used formula Eq. (991 with the two- loop 
running coupling[ Nt = 4 and a s (2 GeV) = 0.31]. We observe that the Sudakov factor provides a 
reduction of around 10% — 25% depending on the choice of A, and changes quite slowly in the given range 
of Q 2 . Therefore we can conclude that the soft spectator scattering contribution may provide quite a 
substantial effect over an extended range of Q 2 if the SCETj FFs fi t % are not too small. 

However, the full next-to-leading evolution includes also single logarithms described by the kernels 
Ua,b- In the ca se of F 2 , the evolution effect from Ub[t \Hh, Hi] in Eq. ( 155 ) is given by the appr oximate 
solution of Eq. (110 1 and does not depend on the gluon momentum fraction r. Using Eq. (110) we can 
write 

F ( 2 S \Q) = ^ e ~ s(Q;Alh ' w) J drU B [t ; Q, fii] /afaQj./ij) 

Jo 



„-S(Q;/j^,M/)rrapp i 



(158) 



Therefore in the ratio F2/F1 the leading and next-to-leading Sudakov logarithms cancel and we obtain 
that this quantity depends only on the SCET/ FFs: 



2 F 2 (S) U* pp [Q, fij] m 2 N dr / 2 (r; Q u ^) 



is) 



U A [Q,m] 



h{Qi,Pi) 



The ratio of the kernels Ua.b in Eq. (159) changes slowly, for instance 



0.93 < U b P ^^ 



< 0.95 for 4 GeV < Q 2 < 16 GeV and A = 400 MeV. 



(159) 



(160) 
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For large Q values, when \Qj\ ~ \/AQ — > oo we expect the asymptotic Q 2 F 2 /Fi — > const as it follows 
from SCET counting rules. It is clear that such asymptotic could be reached only at very large values 
of Q 2 . Therefore it is not surprising that the ratio, measured recently up to Q 2 < 8.5 GeV 2 [7], shows a 
behavior which d rops less fast in Q 2 , when compared with the expected power Q~ 2 . For such values of 
Q 2 the ratio Eq. ( 159[ ) is defined practically only by the ratio of the SCET form factors depending on 
Qi. But the hard-collinear scale in this region is approximately Qi ~ \/AQ ~ 0.9 — 1.3 GeV, i.e. quite 
small in order to expect the asymptotic behavior. 

We obtained that the effect from the Sudakov suppression in the region of moderate space-like Q 2 can 
be estimated as ~ 10% — 25%. However the situation can be different for timelike momenta q 2 > 0. 
In this case, the Sudakov factor after analytical continuation from spacelike to timelike region may 
produce a substantial enhancement. Properties of the timelike processes have been studied in many 
publications (see, for instance, [24l [53] [54]). It is well known that analytic continuation of the Sudakov 
FF to the timelike region produces enhanced it 2 terms. Such corrections were resummed for different 
processes |53H55j . In order to perform such resummation it was suggested to perform the matching 
at a timelike renormalization point — /i 2 [55]. Then the timelike Sudakov factor e~ STL accumulates 
the large ir 2 contributions together with the Sudakov logarithms. Using this recipe we must compute 
Stl = S[—q 2 — ie,—Hh — ie,n] in the timclike region. This can be done with the help of analytical 
continuation of the running coupling which to our accuracy reads [55 , 56 : 

— 2- = 1 - ia[p + —— In 1 - ia{n + 0{a s ), (161) 

where a(fj?) = /3 a s (A*^)/4 ~ 0(1) for moderate values of Q 2 . 

Existing experimental data for the ratio Rm = |Gm(? 2 )|/Gm(Q 2 ) show a considerable enhancement 
of timelike FFs over their spacelike counterparts: Rm — 1-5— 2 over the range Q 2 ~ 10 GeV 2 . Despite 
that the extraction of the absolute value of the time like FF \Gm{q 2 )\ involves considerable assumptions 
about the behavior of the timelike electric FF Ge and probably includes large systematic errors, the 
timelike enhancement is considered as a well established fact. In Ref. 24J it was suggested that "soft 
terms" accompanied by the Sudakov double logarithms could play an important role in a so-called, K- 



factor type enhancement to hadronic FFs in the timelike region. Using the results of Eqs. ( 154 155 ) with 
resummed Sudakov logarithms we can easily estimate such an effect in our approach. 

In order to study the qualitative effect of the SCET/ evolution we consider the ratio of the Dirac 
timelike (TL) and spacelike (SL) FFs. Let us introduce: 

|F!( g 2 )l ^ |e-^||t/J L l LMgliM! ™ 

1 Fi(Q 2 ) " e-Ssr. U A /i(Q/,/i/) ' 1 ' 



where we used qi in order to specify timclike momentum transfer \qi\ ~ yAq in the numerator. We assume 
that the soft spectator scattering mechanism dominates also in timelike kinematics and appropriate 
quantities are related by analytical continuation. At present we do not know the SCET FFs fx %■ However, 



we can study the pQCD evolution effect from the resummed logarithms. In Fig 12 we demonstrate the 
timelike (TL) to spacelike (SL) ratio |e~ STL \/e~ SsL which represents the Sudakov logarithms in the 
ratio Rx- One can see that the obtained ratio |e~ STL |/e~ SsL very weakly depends on the choice of A 
and provides an almost 30% — 40% enhancement effect of the timelike FFs relative to their spacelike 



counterparts. When we combine the Sudakov evolution with the Ua~ kernel of Eq. (100 1 we obtain the 



results shown in Fig 12 (right panel). We see that single logarithms increase the ratio by 5% — 8%. 
But this effect is only a small fraction of the full evolution, i.e. non-Sudakov logarithms cannot provide 
substantial enhancement. Therefore we can conclude that the soft spectator scattering mechanism plays 
an important role in the discussion of the timelike FF's. Sudakov logarithms appearing in this case 
provide an important enhancement in the region of moderate values of timelike momentum transfers 
q 2 . This enhancement is in qualitative agreement with the extracted absolute value |Gj\/|. Moreover, 
taking account of the simple relation of the pQCD evolution in the spacelike and timelike regions, we can 
assume that the enhancement in the timelike region suggests an additional, indirect confirmation of the 
dominance of the soft spectator scattering mechanism in the spacelike region. This might be true if the 
SCET FFs /12 are not modified very strongly after analytical continuation from space like to timelike 
regions. 
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FIG. 12: Left panel: Contribution of the timelike (TL) to spacelike (SL) ratio of Sudakov logarithms je~ STL \/e~ Ssij 
as function of q 2 for different choices of A. Solid and dash-dotted curves correspond to LL approximation with 
A = 400 and 600 MeV respectively; dashed and dotted curves describe next-to-leading approximation. Right 
panel: complete NLL evolution (including kernel Ua) with A = 400MeV. 



VI. CONCLUSIONS 



In the present work, we studied the soft rescattering contribution to the nucleon Dirac and Pauli FFs. 
This work is motivated by phenomenological studies of nucleon FFs suggesting that in the Q 2 range 
5-10 GeV 2 , the nucleon FFs are not yet dominated by a hard scattering mechanism involving three 
active quarks, interacting via hard two-gluon exchange. In the soft rescattering picture studied here, as 
first suggested by Feynman, the highly virtual photon interacts with one active quark, whereas the other 
spectator quarks remain soft. Such a picture is characterized by two large scales : the hard scale Q 2 , 
representing the virtuality of the hard photon probe, and the hard-collinear scale AQ (with A a soft scale 
of order ~ 0.5 GeV), corresponding to the virtuality of the active, so-called hard-collinear quark. 

By way of example, we started our investigation by calculating within perturbation theory the soft 
rescattering contributions to the nucleon FFs. Within such perturbative calculation, the three collinear 
quarks in the initial nucleon wave function are connected to the active hard-collinear quark and two 
remaining soft quarks through (hard-collinear) two-gluon exchange. Analogously, the hard collinear quark 
after the interaction with the hard photon also scatters with the remaining two soft quarks through 
(hard-collinear) two-gluon exchange. For the Dirac FF F\, this analysis overlaps with previous work 
in the literature, whereas for the Pauli FF F 2 it has been performed for the first time here. We have 
demonstrated that such a specific two-loop contribution to the nucleon FFs gives the same scaling behavior 
as the hard region, involving hard two-gluon exchange, i.e., F\ ~ 1/Q 4 , and F 2 ~ 1/Q 6 - Furthermore, 
the perturbative calculation suggests a factorization formula for the FFs in terms of nucleon distribution 
amplitudes, describing how the collinear quarks make up the initial and final nucleon, a hard scattering 
process on the active quark, and a soft correlation function describing the propagation of the remaining 
two soft spectator quarks. 

The specific perturbative calculation demonstrates that a description of the soft rescattering mechanism 
could be carried out in general in two steps. First, one integrates over hard fluctuations (of order Q 2 ), 
leaving only hard-collinear virtualities (of order AQ) and soft virtualities (of order A 2 ). For large enough 
scale Q, such that AQ ^> A 2 , one can then further use perturbation theory and also factorize hard-collinear 
fluctuations leaving at the end only collinear and soft modes, describing the soft QCD dynamics. 

The possibility of such a two-step factorization, with the aim of developing a systematic approach of 
the soft contribution in the case of nucleon form factors, was addressed for the first time in this paper. 
A similar approach has also been considered recently for inclusive cross sections in 57]. The first step 
corresponds with the matching of full QCD onto the soft collinear effective field theory at a factorization 
scale fi 2 = Q 2 , and denoted by SCET/. Technically, we have demonstrated this step by calculating 
the leading-order hard coefficient functions in front of the operators constructed from SCET/ fields, 
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corresponding with the Dirac and Pauli FF structures. These leading-order hard coefficient functions 
involve the emission of hard-collinear transverse gluons, comoving with the active quark. We subsequently 
resummed the large logarithms of order ~ ln<3 2 //i 2 , which appear when evolving the SCET/ operators 
from the hard scale Q 2 down to the scale /i 2 ~ QA. Both for the leading e.m. current operator structure, 
corresponding with the Dirac FF F\, and the subleading operator structure, corresponding with the 
Pauli FF i 7 ^, we solved the renormalization group equations for the corresponding coefficient functions, 
and obtained the NLL solution. This provides a practical check that to NLL accuracy the first-step 
factorization (so-called SCET/ factorization) for both F\ and F2 indeed holds. 

We next discussed the further matching of the SCET/ theory to the effective theory involving only 
collinear and soft particles (so-called SCET//), defined at a factorization scale \j? u = QA. As a first 
step to arrive at such a full factorization formula for the soft rescattering contribution, we analyzed 
in this work the leading terms in the effective theory. The factorization formula involves two so-called 
jet functions, describing the amplitude for the transition of three collinear quarks into a hard-collinear 
(active) quark and two soft quarks; a soft correlation function describing the soft rescattering of the two 
soft spectator quarks in the background soft-gluon fields emitted by the hard-collinear (active) quark; and 
the two nucleon distribution amplitudes, describing how the three initial and final collinear quarks make 
up the nucleons. The jet functions can be computed performing the matching from SCET/ operators 
onto SCET// at the factorization scale p? u = QA. Also here large logarithms ~ lnAQ/^jj arise, when 
we evolve the factorization scale fi 2 j down to value of order A 2 . They can be resummed again using RG 
equations. We leave this consideration to a future work. 

For the Pauli FF F% we also discussed that an analysis is more involved as there may be a double 
counting between the hard and soft rescattering mechanisms. Furthermore the matching from SCET/ 
onto SCET// does not yield a well defined expression for the Pauli FF, due to end-point singularities, 
which calls for a more refined treatment for F 2 in a future publication. 

The SCET/ factorization formulas allowed us already to discuss some phenomenological consequences 
in this work. For the soft rescattering contribution to the Q^Fn/Fi ratio, we found that the ratio of the 
next-to-leading order evolution kernels changes only by a few percent in the range Q 2 ~ 4 — 16 GeV 2 , 
and is mainly dominated by SCET/ FFs defined at a corresponding scale Q 2 ~ AQ ~ 0.8 — 1.6 GeV 2 . 
Such scale is quite small to expect the asymptotic constant behavior. The experimental data for the 
Q FijF\ ratio in this Q 2 range indeed show rising behavior, in agreement with the above analysis. A 
second phenomenological consequence of our framework was discussed for the ratio of the spacelike to 
timelike FF F\. We showed that the resummed Sudakov logarithms provide a 30% - 40 % enhancement 
to this ratio in the range of momentum transfers q 2 around 10 GeV 2 . This enhancement is in qualitative 
agreement with the empirical extracted ratio for the absolute value of the dominant FF Gm in the 
timelike as compared to the spacelike region. A more detailed phenomenological analysis requires us to 
parametrize the SCET/ FFs, which is equivalent to using the SCET// factorization formula, and express 
them in terms of DAs, jet functions, and a two-quark soft correlation function, as outlined in this work. 
Such an analysis remains a challenge for a future work. 

Appendix A. Leading order coefficient functions 



Here we discuss in detail calculation of leading-order hard coefficient functions. First, Eq. (87 1 can be 
rewritten in compact form in momentum space, 

g(0)7^(0) = tr[C"(Q,/i) O*(0)] + [ dr tr [ C£(t,Q,h) Ol{r) + C£(r, Q, //) 0«(r) ] + ... (163) 

Jo 

where we used translation invariance and defined the momentum space coefficients as 

Cf*(Q, At) = / dhdh G>(s!,i 2 ) e i{P'-n)8*-i{P-n)s^ (164) 



C^Q^) = I dhdhdh e~ l(p - n >i e l[p ' ^ «»+*('"•«) <5£(s!, s 2 , s 3 ), ( 165 ) 
CZ(t,Q,v) = J dMhdh e« p '- n ^ s 2 , s 3 ), (166) 
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Here P and P' denote the total hard-collinear momentum of the external state for each jet, and (P ■ n) = 
{P' ■ n) — Q is the large component of each momentum. The variable r is the fraction of large momentum 
component (P' • n) [{P ■ n)] carried by the hard-collinear gluon A'j_ (A±), f = 1 — r. The objects 0|(r) 
and 01(t) denote the Fourier transformed SCET operators: 



OUr) ~ 



ds 3 
2ir 



e -« 3 (P'-«)ro| (0)0)S3)) q , (t) = 



ds 3 Js 3 (P-n)r 

2tt 



O g n (0, 0,8 3 ). 



(167) 



The tree-level coefficient functions in momentum space can be obtained from an analysis of the matrix 
elements in QCD and SCET. In order to compute C^(Q,fj.) defined in Eq. (1631 consider the matrix 
element of the e.m. current between collinear quark states: 



(p'\mrq(o)\p) = tfMCiQ,n) o 9 (o)]b> 



(168) 



The subleading term in Eq. ( 163 ) does not contribute in this case. Then for the matrix element at LO 
we obtain: 



IhsoiEq. ((T68J): (p / |?(0)7^(0)b> Ij o=«V" = f7^ + O(l/Q), 



(169) 



rfeofEq. p§): (j/ |tr [<7"(Q, ft) O q (0)}\p) = tr [C"*(Q,/x) £' ® £] , 



(170) 



where without subscript ft.c denote large components of Dirac spinors (|30j) and ( |3l[ ). Comparison 
(fl69l) and (fl70l yields 



(171) 



where F',F describe quark color indices. 

In order to compute the subleading coefficient functions one has to consider the matrix element with 
the quark-gluon external state. We consider an outgoing gluon with hard-collinear momentum q' collinear 
to p' and for simplicity we neglect the transverse components of the outgoing momenta. Then we can 
compute the leading-order contribution to C^(r, Q,/i). 

We start by considering the QCD calculation. The corresponding diagrams are shown in FigfTB") For 





q , a, v 



FIG. 13: QCD tree diagrams required for matching onto the subleading 3-particle operator 



the first graph we have 



2(p' + q')-n 2(p'+q')-n 
± (p' + q') ■ n\ 



(172) 



For clarity, we write for the external gluon line with momentum q' instead of polarization s*(q'). The 
second diagram: 



d 2 = (i ff )u(p')7 /t ^— S^'«(p) * da) 

[p-q'y 



(ig) 41 



tf> l 7/1 l 

2-{q' -n) + 2 • n) 



(p ■ n) 



(173) 
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Therefore the sum reads 



Di + D 2 = £' 



(~9) 41 



(p ■ n) (p' + q') ■ n 



(174) 



which involves the transverse gluon field and the longitudinal projection of the e.m. current, as required. 
It is easy to see that the obtained kinematical structure is e.m. gauge invariant. This term must be 
compared with the SCET matrix element, 



dr tr[ C£(t,Q,h) 0|(t)] 



ith 



(p', q '\OUr)\ P )^eA^ S 



P )= I dr tr[C£(T,Q,A0 (j>',<f\OUr)\p)}, 



(p' + q') ■ n 



Substituting this into Eq. (175) and comparing this with the QCD result of Eq. (174 1 we obtain 

1 



C£(t,Q,m) 



Q 



(n" + n") 1 5 



F'F, 



(175) 



(176) 



(177) 



where the symbol 1 denotes the unity operator in Dirac space. Notice that the obtained coefficient 
function does not depend on the momentum fraction r at the LO level. A calculation of the second term 
with can be done in an analogous way. The result can also be obtained without explicit calculations 
by invoking time reversal invariance which demands the result to be symmetric under n, ■<->• n. 



Appendix B. Correspondence between QCD and SCET calculations 

In order to illustrate the correspondence of SCET with QCD we here perfor m th e calculation of the 
jet functions discussed in Sec. [n] We start from the calculation of J' defined in ( 122 ). In order to have a 
direct correspondence with the expression ( 48 1 



v r i 11 

J (a) [Vi^i] = 



Q 3 yivl ( w i + w 2 ) 2 (-wi) 

we consider the appropriate subprocess 

&c(p' -ki - k 2 ) + q s {k 1 )q s (k 2 ) S °4- Tl £ c (Vlrf) tcivip') £c(V3P')- 
described by the matrix element 



J{ a) QaMQaM = (yip'^p'mp'lT^W'iO),^,^ 



(o) 



cr ) io) 



(178) 
(179) 



The soft quark fields q s are considered as external. In order to reproduce the expression in Eq. ( 48 ) we 



need the diagram shown in Fig 14 



From the Lagrangians and C^ 1 one can easily define the Feynman rules. They were already 
presented in the literature (see, e.g., Refs. [571 EH HZ] ) • For the convenience of the reader we reproduce 
the relevant vertices here. Taking account of only the required leading-order terms one obtains 



^(1) 



fi,a 



[i, a 



* h 



ig[T*\ 



AB 



(180) 



jx, a 

4lhc 




v, b 
4 2 hc 



-19 



[T a T b ] AB hlU 



(Q2 ■ n) 



[T a T b ] AB hlU 



(qi ■ n) 



(181) 
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y,p- 
y 2 p'- 
y 3 p' 



FIG. 14: One of the diagrams describing subprocess 
fields, and outgoing quarks are collinear. 



1781) in SCET/. The soft fields are considered as external 



Assuming the same choice of momenta as in Fig. [3j we obtain the following analytical expression 

1 1 1 



" '_L rpbrpa 



1 



(Pig ■ n) 
1 11 



T b T a \ q{k 2 W-r {T a n»} , , 

' 4 (p 3g -n)pj g pl g 



(k+ + k+Y (- fc+) Q 3 yivl 

d^i,2 J[ a ) [Vi,Ui] q ai (ki) q a2 (k 2 ) 5{u) X - k^)5(uj 2 - 



(182) 



with the same J',-, as in Eq. (48 1. Consider the color structure which we ignored in the calculation in 
Sec.lnJ Projecting the color indices of the outgoing collinear quarks onto the colorless nucleon, we obtain : 



s i'j'k' 

3! 



[ T %i® [T b T a ].,.®[T a [ 



k'k 



Jjk 
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(183) 



The resulting antisymmetrical tensor e % ^ k is then contracted with the color indices of the soft fields 
yielding the soft operator in Eq. (131). 

Consider now the helicity flip FF F 2 . Again, from the QCD calculation we obtained the result of 



Eq. (681: 



J [a) [Vi,Ui] 



2/12/3 (W2 + w 3 )w| 



(184) 



In this case, we have to consider the presence of the small component 773. In SCET this field is eliminated 
by the equation of motion, yielding 



In order to have a connection with the QCD result of Eq. (68), we must substitute 

ff c = &D ±c (infi c )-^ ~ ? c A c± (in ■ V)" 1 * 



(185) 



(186) 



From this expression one can see that such a state consists of a collinear quark and a transverse gluon, 
as is shown in Fig. [9J The corresponding vertex is generated by the Lagrangian C^j , includes 2 gluons, 
and can be associated with the following combination: 



£ hc W±Wq ~ cj c gA' ±c 



(fid) 



ign ■ A hc 



ffj! ±e (n- d)~ 1! ~ 



A hc-^- q{k 3 )- 



(187) 
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Therefore, we obtain for the diagram in Fig 14 (again ignoring color structures) : 

jjll/l '. 



Q 2 ym (k+ + k+) (k+) 2 
= / rfwi :2 J„ [yi,Ui] q ai (h) q a2 {k2) 5(u>i - k^)6(uj 2 - fcj), 



(188) 



with the same as in Eq. (184). 

Using these two examples, we demonstrated that the SCET correctly reproduces the tree-level hard- 
collinear subprocesses computed before in QCD. 
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